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mercury itself can then be dried with bibulous paper, and after 
being gently heated becomes pure and bright. 

20 . I will only add,' that the modifications which these experi¬ 
ments suggest in the construction of the ? common portable artificial 
horizon are very obvious. With regard to the ingenious arrange¬ 
ments which, after? the completion of these experiments, Colonel 
Sabine informed me had been adopted some years ago by M. 
Schumacher for jthe artificial horizon, I found the formdf the vessel, 
viz. the segment of a sphere, peculiarly objectionable, on account 
of the depth and the quantity of mercury employed: moreover, 
there is no possibility of easily cleaning the surface, or of removing 
the obnoxious film of amalgam so often alluded to. 


On the Determination of the Orbit which a Satellite describes 
around its Planet By M. Sawitsch, Professor of Astronomy 
in the University of St. Petersburg; 

When there exist a great number of observations of a satellite, 
we may employ various combinations for finding with facility the 
elements of the orbit which the revolving body describes around its 
primary. It appears, however, that no method has yet been as¬ 
signed for determining these elements from the. smallest possible 
number of observations. The general theory which we owe to the 
illustrious Gauss, relative to the calculation of planetary orbits, 
undoubtedly applies also to. the problem under consideration. 
This theory first leads to an approximate solution; after which it 
puts us in the way of obtaining such a degree of accuracy as the 
excellence of the observations demands; The provisional results, 
however, necessarily imply that the intervals between the times of 
observation are small, relatively to the complete revolution of the 
celestial body in its orbit. Such a supposition cannot be admitted 
when the question refers to a satellite, for in that case a single day 
may form a considerable part of the entire period of revolution. 
It is proposed, according to the following method, to determine the 
elements from three observations, the interval between the first and 
third being supposed not to exceed a few days.: 

Let f, t 1 , t”, denote the times of the first, second, and third ob¬ 
servations; A, A', A", the observed distances of the satellite from 
the planet; and Q, Q', Q", the angles of position referred to the 
circle of the geocentric latitude of the planet. Moreover, let 
7, l', Z", b y b’ y b” t denote the geocentric longitudes and latitudes of 
the planet: and R, R', R", its corresponding linear distances from 
the earth. 

All these quantities are supposed to be given either by ob¬ 
servation or calculation; but the following are to be determined, 
viz.: — 

R (i + §), R' (i + £'), R" (i +<T), the linear distances of the sa¬ 
tellite from the earth, corresponding to the times t, t , t' r . 
r, the period of the. revolution of the satellite. 
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the angle under which the radius of the orbit of the satellite 
would appear at the mean distance of the planet from the 
sun. 

t, the inclination of the orbit of the satellite to a plane passing 
through the centre of the planet parallel to the ecliptic, 
the longitude of the ascending node of the orbit measured 
upon the same plane. 

f, the longitude of the epoch. 

Let the motion of the satellite, as well as that of its primary, be 
referred to a system of rectangular co-ordinates, the origin of 
which is supposed to be at the centre of the earth. Let us suppose 
that the axes of X and Y are in the plane of the ecliptic; also let 
the positive axis of X be directed to the longitude Z 7 , the positive 
axis of Y to the longitude 90°+ V, and the positive axis of Z to the 
north pole of the ecliptic. 

Let X, Y, Z, X 7 , Y r , Z',X", Y", Z", be the co-ordinates of the 
planet, and X -j-x, &c., X 7 -}- x', &c., Z"-f z u , &c., those of the sa¬ 
tellite corresponding to the times t , t', t”. Let d l, d l', dl ", repre¬ 
sent the differences of geocentric longitude of the planet and sa¬ 
tellite, and db, db', db", the differences of latitude. Hence we 
have" • 

X = R cos b cos (l — l), X + oe — R (1 + 2 ) cos (b + d b) cos (l —l'+ d l ) r 
Y = R cos b sin (l — V), Y + y = R (1 + S) cos (b + db) sin (l — l'+dl) 
Z =» R . sin b f Z + z =■• R (1 + V) sin (b + d b ). 


Similar expressions may be obtained for the other co-ordinates 
of the planet and satellite. 

We have also manifestly 


dl 


A. sin Q, 
cos b ’ 


dl' 


A'. sin Q/ 
cos 6' ; 


dl" 


A ,r . sin Q" 
cos b" 5 


db = acosQ; db'=* A'cosQ'; db H ~ A "cos Q,". 


Let us now transport the system of co-ordinates parallel to 
itself from the centre of the earth to the centre of the planet. The 
co-ordinates of the satellite relative to the centre of the planet will 
then be represented by x, y, z, x', y', z, x”, y", z . Remarking that 
R*, R'*, R" y, sin dl, sin db, sin dl, sin db', sin dl', sin dbd 
may be considered as small quantities of the first order, the squares 
and products of which may be neglected; and bearing in mind 
further that l l, l" — l, do not amount, during an interval of 
several days, to more than a few minutes of space, we have without 
sensible error, 

x = R . S cos b — R . sin b . cos Q,. A . sin 1" *] 

y — R . sin Q . A . sin 1" l (A) 

z = R . $ sin b + R cos b . cos Q . A . sin 1" J 


Similar expressions may be obtained for xy', z', x", y", z . 
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Let n denote the area of the triangle formed by the radii vec- 
tores of the satellite corresponding to the second and third ob¬ 
servations, with the cord which unites their extremities. Let n ' 9 n" 
represent the areas of the similar triangles corresponding to the first 
and third, and to the second and first observations. Since the satel¬ 
lite revolves in a plane passing through the centre of the planet, 
which coincides with the origin of the co-ordinates x 9 y , z , the fol¬ 
lowing equations are deducible:— 

+ »o (i) 

n y — n! y* + n” y” *= o (a) 

nz--n'z'+ri f z".— o (3) 

Substituting in (2) the values of y 9 y , y ", as already determined, 
we have 

B • A ♦ sin Q l£. r "- a "- sin ( \ 

n ~ R'.a' . sinQ' * n' R' . A' . sin Of 

Multiplying (1) by sin b\ and (3) by cos b\ and subtracting 
the one result from the other, bearing in mind at the same time, 
that as the geocentric latitudes of the principal planets are gene¬ 
rally small, their alterations, b — b\ b" — b', in the space of a few 
days are inconsiderable, we have 

n' R.a.cosQ n" R". a", cos Q!' 
n ~R'. a'.cosQ' + n ' R'. A' . cos Q' ^ 


If we.suppose the satellite to revolve in a circular orbit whose 
radius is r, we may employ the equations («) and (fi) in finding r 
the period of revolution. 

Let v”, v 9 v, be the angles described by the radius of the satellite 
during the intervals of time t, t” — t, We have then 


n 


-i- r 2 sin v, n' — — r 2 sin «/, n" 
2 2 


r « . if 

■— r 2 sin v ' 
2 


rtf sin v f n" sin v" 

— = . —l* —= —- 

n sin v n sin v 






360° 





(f - *) 


It is clear, then, that if an approximate value of r was known, we 
might easily find v, v f 9 v " 9 and hence also the relations —, —• Now 

it is always possible to obtain such a value of t by successive 
trials. Assuming r = 3 or 6 or 10, &c. days, we calculate the cor¬ 
responding values of s ! n - ~, Sm V - , and then we shall see which 
r ® sin V sm v 

of the hypotheses satisfies best the equality of -^7-— and — as 

determined by (*) or (/3). The corresponding value of r will express 
the approximate value of the period of revolution. In order to avoid 
useless calculation, it will be desirable to take into consideration 
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the three angles of position Q, Q', Q", which will serve to show if 
the interval of time t'' — t between the first and third observations 
approaches nearly or not to the time of a complete revolution. 
Knowing r approximately, and repeating the calculation fox two 
hypotheses of that element, differing slightly from each other, we 
may obtain a more accurate result, and in this manner we may 
approach nearer and nearer to the true value. When the angles 
Q, Q', Q", are very small, or very nearly equal to i8o°, it will be 
advantageous to employ equation (/ 3 ) instead of equation (*). 

The following example, based upon three of M. Otto Struve’s 
observations of Neptune’s satellite, will serve to illustrate the prac¬ 
tical utility of the method:— 


Date. 

M. T. Pulkowa. 

Angular 

Distances. 

Angle of Position 
with respect tb the 
Circle of Latitude. 

1847. Sept. 25 

IO h 2 I m = t 

15^92 = A 

228° 45' = Q 

28 

I 2 b 9 m = t' 

17"* 30 = a ' 

57 ° 4 *' = Q' 

30 

9 h 36“ = t" 

6 

II 

l> 

164° 51' = Q" 

Hence 

t = 3 d *o 75 , t" - f = 

i d *894, t" - / 

= 4 d * 9 ^ 


Moreover, the ephemeris of Neptune gives 


log R = 1*46450, log R' = 1*46505, log R" = 1*46540 

It will be seen at once that the time of a complete revolution 
exceeds five days, for in 4*969 days the angle of position (Q") has 
not nearly resumed its original value (Q). We may then com¬ 
mence by supposing r = 5*5 and r = 6. The former of these 

hypotheses gives for ■ ^ ■ a numerical value by far too small in 

n * 

comparison with its equivalent — calculated from the formula (#)♦ 

• si n ty* 

Taking r = 6, we obtain for —:-a value too great. If we now 

° ’ sin v * 

assume r = 5*7 and r = 5*8 we obtain 


. n' ■ sin v f , , N 

log-log —:- = + 0*00643 (r = 5*7) 

n sin® 



log 


sin v 
sin v 


0*02390 (<r 


5 * 8 ) 


Thus the one hypothesis is 0*00643 in excess, and the other, 
0*03390 in defect, and hence by proportional parts we find 
T = 5*72. Repeating the calculation with r = 5*72, we shall find 
that this hypothesis-is erroneous only to the extent of + 0*00113, 
and we obtain definitively r = 5*724 days. By a discussion of the 
totality of eighteen observations made in the years 1847-49, M. A. 
Struve found t = 5*876 days, a result which accords sufficiently 
well with that deduced by an approximative process from three 
isolated observations. 
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When once we know the value of r we may easily determine the 
radius of the orbit of the satellite. 

Let r denote the linear value of this radius; A the semi-major 
axis of the orbit of the planet around which the satellite revolves ; 
v" the angle included between the two radii whose extremities are 
determined by the co-ordinates x, y,z, x', z\ We have then 

' - CGS V ,f == 

substituting for x, y,.z, x 9 y’, z! their values (A); bearing in mind 
also that g = A ^ and that r 4 = x % -j- y* *+• z % = x lZ -f- y fl + z n , 
we finally obtain 


.« =/fK*±v/fK 4 - M 1 

in which 


M ! = (x A ) (t - 4 )' sin2 ( Q,_ Q ) 

sin 2 «" . . ; ... : 

R* = (t- 4 ) s+ (t■ A, )* _ 2 (x• A ) (x- A ') cos(Q - cos 

■ sin 2 v" 


It is easy to select the proper sign of the quantity under the 
radical in the expression for g, since the latter quantity is essentially 

R A R / A / R // A /; J 

greater than —f—, —j —, or —j—. We may determine g in a 


similar way by combining the third observation with the first or 
second. The most favourable combination is that in which sin v 
has the greatest value. 

In the above example we have, 


Q«= 228°45 / Q / =57°42' 

v" ~ 193 0 22'; -S. . a = Is' , *445 ; ~ . a' = i6"-8o6 

l K 2 = 206*3 ; M 2 = 30624; 109*2 

e 2 = 3i5*5 5 ? = i7"*76a 

M. A. Struve has deduced g = i7"*969 from the totality of 18 
observations. 

We have now to calculate the planetocentric longitudes and 
latitudes of the satellite, A, A', a", 6 , 6 ", with respect to a plane 

passing through the centre of the planet parallel to the ecliptic. 
The expressions of the three planetocentric co-ordinates x 7 y,z, in 
which the direction of the axis of a? is determined by the longitude % 
lead to tire following equations 
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R 

g . cos 0 . cos (x - /) = | cos 6 — sin b -r- . A , cos Q («) 

A. 

R 

e . cos 0 . sin (x — 7) = — . A . sin Gt ' ( b) 

R 

£ sin 0 ~ | sin b + cos b . . A . cos Q (c) 


where i = ± ^e“-(f- • *)* 

Dividing (6) by (a) we find an expression for tan (a — Z), and 
as l or the geocentric longitude of the planet for the time t is 
known, we thus obtain the value of A. The equation (c) serves to 
calculate 0, the planetocentric latitude of the satellite. 

If we suppose the direction of the axis of x to be determined by 
the longitude l ', and then by the longitude l", we shall obtain equa¬ 
tions analogous to (a), ( b ), (c), which will serve to calculate A 7 , 0 ', 
and a", 0", by means of the data relative to the second and third 
observations. 

When we have determined A, 0 , a 7 , 0 ', the following formulee 
will serve to calculate the inclination i, and the longitude SI of the 
ascending node of the orbit of the satellite upon a plane passing 
through the centre of the planet parallel to the ecliptic. 

sin . tan i — \ (tan # + tan 0) sec \ (x 7 — x) 

cos . tan i = | (tan 0 7 — tan 0) cosec \ (x 7 — x) 

We might obtain i and ^ in a similar manner by combining 
A 77 , 0 ", with A 7 , 0 ', or with A, 0 . 

In our example, 

o / o / 

by combining x', 0', with x, 0, we have SI, = 2.64 46 ; « = 37 35 

— — ; X 77 , 0 " — X 7 , 0 ' — & = 296 36; i = 30 28 

The difference of these two determinations, although sufficiently 
great, does not depend upon calculation. It arises from very small 
errors to which the best observations are liable, and which may 
exercise a very sensible influence according to the positions of the 
satellite: with respect to the nodes at the times of observation. We 
may remark that the combination of the values A, 0 , A 7 , 0 ', is not 
favourable for determining Si, since the satellite was in those 
positions nearly —90° and +90° from the node. The combination 
of the second observation with the third is the most favourable for 
thispurpo.se. As regards i both combinations are almost equally 
good. We admit, then, 

•SI = 296°, i =* 34 0 

With these results, the equations 

tan u = tan . (x — Si) sec i ; A— u 4 Si 
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will give the arguments of latitude u y u\ u \ and the planetocentric 
longitudes A, A 7 , A' 7 , of the satellite in its own orbit. 

In our example, A = 201 0 34', A' = 34 0 o'. A" == 148° 16'. 
Designating by p the mean daily motion of the satellite, we have— 

A n — A ' : ^ 

ft = -y, -— = 6i°* 72 by the interval of nearly 5 days. 

A'-A £ £ ' ' : 

f 4, = f ~ t ~ “ 62 6 4 — — 3 — 

A" - A' ■ ' 

* = T—7 ~ 6o ’35 - - - • - 

Whence the probable value of the mean daily motion, 

(a =. 61 0, 72, 

and the period of revolution, 

«r ss 1^2- = 5*832 days = 5 d I9 h io m . 

With p = 6i°* 72 we find the planetocentric longitude of the 
satellite upon its orbit, in 1847, September 28, at I2 h tf 1 mean 


time to be — 

By the first observation and the interval of 3*073 d .. 31 0 14/ 

second — • - —— o . 34 o 

third — — 1*894 . 31 26 

Mean value . 32 13 


The corresponding argument of latitude is 96° 13'. Subtract¬ 
ing the motion in i roi days, we have 136° 45' for the argument 
of latitude onv the n-th of September, at n h 55 111 mean time. 
M. A. Struve determines the same argument, by means of eighteen 
observations, to be equal to 128° 26'. We are, therefore, in error 
only to the extent of 8° 20'. 

As the inclination of the orbit of Neptune to the ecliptic is 
= i p 47V and as the longitude of the ascending node of the orbit 
is = 130° 7', we hence deduce the longitude of the ascending node 
of the satellite upon the orbit of Neptune to be = 296° 37', and 
the inclination to be = 35 0 44'. We have now to explain the modes 
of verifying what has been already found, and to assign results of 
greater accuracy. 

Since the trustworthiness of the preceding calculations depends 
upon the precise knowledge of r the period of revolution, which 
has only been determined approximately in the first instance, it will 
be desirable to ascertain whether this value accords sufficiently 

well with the value r = The latter value ought to be in 

ft* 

general the more accurate of the two, and if it differs considerably 
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from r as determined by means of («) or (&), it will be necessary to 
repeat the calculations. 

Assuming p = 6i°*J2, we have r = 5*832 days. By the first 
calculations we obtained r = 5*724. The difference does not 
amount to more than o* 108 days. 

Let us now inquire into the degree of approximation of the; 
angle g as already determined. With y, = 6i°*72 we have-— 


' = 189° 40' } Q' - Q = 188 0 58?; 2^ -1 5"*445 5 


R' A 7 
A 


= i6"*8o6 ; 


whence £ = 20" *40 


v = 116 0 54'; Q" — Qf = 107° 8'; = 4*929 ; 


R' A' 
A 


i6"‘8o6; 


whence £ = 18*27 • 

r« '^/ / R A 

*' - 3 o6 ° 4 2 '; Q'- Q =- 2 ’9 6 ° 54 '; —= 4 ' 9 2 9 » = J 5'445 ? 

whence g = 18*02 


As the. errors A, A', Q, QY produce upon the determination of g 
an influence so much greater, as sin v” is less, we may take—- 


20"*40 sin v " 4 - i8"*27 sin v 4 * i8"*02 sin v' __ j$f'. 

sin t>" 4 sin v 4- sin V 1 ^ 

for the probable value of g, as deducible from three observations. 
A new calculation of SI and i would diminish only in a small degree 
the uncertainty of these elements. 

It remains for us to say a few words upon the mode of obtaining 
the elliptic elements of the satellite. Although four observations 
suffice theoretically for that purpose, still, considering the small 
value of the excentriqity and the sensible influence of errors of ob¬ 
servation, it will be impossible to attain the object by means of these 
data alone. Before undertaking the calculation of the elliptic ele¬ 
ments it is desirable to ascertain the circular elements, and that to 
a certain degree of precision, so as to treat their corrections as 
small quantities of the same order of magnitude as the eccentri¬ 
cities. It is necessary, besides, to have a long series of good obser¬ 
vations, in order to deduce from them the. elliptic elements of the 
orbit in a trustworthy manner. 

The first thing to be done is to ascertain the degree of accuracy 
with which the observations are satisfied by the circular elements. 
For this purpose, it will be necessary to calculate with these ele¬ 
ments the angle of position Q 0 and the apparent angular distance 
A 0 for the time of each observation. Representing by A, Q, the 
values of these quantities as deducible from the corresponding 
observations, we then form the- differential quantities A — A 0 , 
Q — Q q , which will serve to exhibit the errors of the existing hypo¬ 
thesis respecting the orbit of the satellites Now if the orbit be 
supposed to be in reality an ellipse, A — A 0 and Q — Q 0 may 
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generally be expressed in terms of the corrections of the circular 
elements and the excentricity and longitude of the perihelion. Let 
e be the excentricity and ?r the longitude of the peri-planet; then, 
putting e cos 5r = 4 ' and e sin w = n, we may express A — A 0 and 
Q — Q o in terms of seven unknown small quantities, viz.: the cor¬ 
rections of the five elements already involved in the circular hypo¬ 
thesis and the quantities ^ an ^ An equation of condition 
involving these quantities may be formed for every observation, 
whether of distance or position; and from the totality of such 
equations, the values of the unknown quantities may be deduced 
by the application of the method of least squares. 


Extracts of a Letter from M . Otto Struve to the Astronomer 
Royal , dated January , 1853. 


“ As to what regards Vitlarceau’s remarks on the differences in 
right ascension of feeble comets, as observed by different instru¬ 
ments and methods, I am quite of opinion that similar differences 
must exist, and that for different reasons, but principally produced 
by the different optical power of the telescopes. The nature of 
these differences is, perhaps, not yet sufficiently explained for all 
different methods of observation ; but, generally speaking, the most 
powerful telescopes must furnish the most exact results; and for 
this reason, I think, Villarceau has done wrong to attribute a 


higher precision to the Konigsberg observations than to those made 
with the Northumberland telescope, especially as Wichmann him¬ 
self thinks the accuracy of his observations of D’Arfest’s comet 
has seriously been injured by the faintness of the object, while 
Challis does not mention at all that these observations have given 
him particular trouble.” 

“ Last month the University of Dorpat celebrated the fiftieth 
anniversary of its renewed existence. On this' occasion I wrote a 
dissertation, entitled Narratio de Pardl taxi SteUa o& Lyra. I n 
this tract I have given, besides some general remarks on the deter¬ 
mination of parallaxes, the result of fifteen months’ observation on 
the parallax of * Lyra. The comparison star used by me is 
identical with that of mv father; but the result does not so well 


agree with the former, as I thought it would do. My observations, 
which, as to accuracy, appear to be considerably superior to those 
of my father, make the parallax by the 


Angles of position = +0*145 — 0*016 - 
By the distances *= +0*140 ± 0-025 " 

M In combining these results with the former values obtained 
by W. Struve and Peters, the parallax of et Lyra is found in the 
mean ~ 154.9 di o"*oi20. ; j 

“ I have indicated in the same tract several reasons, which will 
account in part for the larger parallax found by my father, but the 
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